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ABSTRACT 

We have measured the dayside spectrum of HD 189733b between 1.5 and 
2.5 fim using the NICMOS instrument on the Hubble Space Telescope. The 
emergent spectrum contains significant modulation, which we attribute to the 
presence of molecular bands seen in absorption. We find that water (H2O), 
carbon monoxide (CO), and carbon dioxide (CO2) are needed to explain the 
observations, and we are able to estimate the mixing ratios for these molecules. 
We also find temperature decreases with altitude in the ~ 0.01 < P < ~ 1 bar 
region of the dayside near-infrared photosphere and set an upper limit to the 
dayside abundance of methane (CH4) at these pressures. 

Subject headings: planetary systems — techniques: spectroscopic 



1. Introduction 

HD 189733b is a transiti ng hot- Jupiter planet in a 2.2-day orbit around a K2V stellar 



primary (IBouchy et al. II2005I ). Due to the relatively large depth of the eclipse (~ 2.5% at K 
band) and the bright stellar primary (Kmag = 5.5), this system was immediately recognized 
as an important target for atmospheric characterization observations, and its emission has 
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been studied extensively durning secondary eclipse (IDeming et al 



2007: Grillmair et al. 2007: Knutson et al. 2008: Charbonneau et a. 



2005 



photometric obs ervations revealed the pres ence of H2O (ITinetti et al. 



Knutson et al. 



20081). Multicolor 



2007af ) and the likely 



presence of CO ( Charbonneau et al. 2008[). while o ptical transmission spectroscopy suggests 
scattering by high altitude haze (jPont et al. II2008I ). In ad dition, extensive theor et ical work 



has a ls o been done on the a tr nosphere of this planet by iFortney et al. I (120061 ): iBarman 



(j2008|); Burrows et al. I (120081 ): IShowman et al. I (l2008f ) and others. 



Recently, the NICMOS camera on the Hubble Space Telescope (HST) was used during 
the primary eclipse of HD 189733b to obtain a near-IR transmission spectrum of the planet's 
atmosphere; these results showed the presence of H2O and CH4 (Swain et al. 2008; hereafter 
SVT08). The near-IR transmission spectrum of HD 189733b probes the upper (P ~ 10~^ 
bar estimated from our models) regions of the of the atmosphere at the terminator. In this 
paper, we report the results of HST observations of the dayside spectrum (derived from 
secondary eclipse measurements). At near infrared wavelengths, the dayside portion of the 
atmosphere of HD 189433b is probed in deeper (P ~ 0.1 bar estimated from our models) 
regions. 



2. Observations 



We observed HD 189733 with the Hubble Space Telescope for five contiguous spacecraft 
orbits, using NICMOS in imaging-spectroscopy mode with the G206 grism (wavelength cov- 
erage 1.4-2.5 fim). Observations began on April 29, 2007 at UT 23:47:58 and ended at UT 
06:51:52 on the next day. The first two orbits (Oi & O2) observed the target pre-ingress, the 
third orbit (O3) was phased to capture the occultation, while the fourth and fifth orbits (O4 
& O5) were post-egress (see figure 1). A total of 638 usable snapshot spectra were acquired 
during the five orbits. The third orbit (the oc c ultation) containe d 130 spectra. Using the 
best available ephemerides (jWinn et al. I l2007l : jPont et al. 1 120081 ) we determine that fuU- 
occultation starts ~ 13 exposures into O3, and lasts for the remainder of this orbit. The 
effective exposure time for each spectrum was T = 1.624 s; the overall instrumental configu- 
ration, including the location of the spectrum on the focal plane array, was identical to that 
in SVT08. A few spectral calibration exposures (in Oi) were acquired using a narrow-band 
filter. 



Because the NIC-3 camera is severely undersampled and because of large gain drop-off 
at the edges in the NICMOS detector pixels, the instrument was configured in DEFOCUS 
mode (FWHM ~ 5 pixels). This defocus reduces the level of photometric fluctuations due to 
both pointing jitter (random) and beam wander (systematic errors). While the defocus also 
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helps to minimize the bright-source overheads for this extremely bright target, it nevertheless 
limits the effective spectral resolution to i? ~ 40. 



2.1. Data Analysis 

A complete description of the data-analysis methods is provided in the supplementary 
information to SVT08. Herein, we discuss only departures and addenda to these analysis 
methods. 

First, each image is reduced to a one-dimensional spectrum covering 1.5-2.5 /im in spec- 
tral range. This includes combining data from the grating's first and second orders, which 
are simultaneously imaged onto the detector array. Inclusion of the second order improves 
the SNR in the blue because the second order, which is partially imaged onto the detector, 
contains more detected photon-fiux at the shortest wavelengths (shortwards of 1.7 /im, long- 
wards of which the efficiency of the 2"'^ order falls rapidly) as a consequence of the peculiar 
blaze of the parent grating. Prior to combination, wavelength solutions for both orders were 
obtained using the calib ration exposures and the wavelength-scale coefficients in Nicmoslook 



( iFreudling et al. 1119971 ). Merger of the spectral orders takes advantage of the near integer 
ratio, ~ 2, of the spectral resolutions of the two orders. The orders are coadded around a 
fiducial wavelength, and any smearing is captured by assigning the coadded channel a wave- 
length equivalent to the weighted average of wavelengths of component channels. Prior to 
coadding, the 2"''^ order spectra were multiplied by a broad Tukey window so as to suppress 
the effects of (1) the sharp band-edge at 1.5 fim and (2) the detector array edge at ~ 1.85 /im 
on the final spectrophotometry. Lastly, coadded spectra are rebinned to wavelength chan- 
nels at the binwidth of the DEFOCUS full-width-at-half-maximum, guaranteeing adjacent 
channels have independent spectral content. As the state variables are channel independent, 
the two orders can be coadded prior to decorrelation (see below). The secondary eclipse is 
observable in the raw light-curves (see Fig. 1) and especially so in spectral subbands with 
high contrast {e.g. around 2.3 fim; see Fig. 3). For example, if we divide the "red broad- 
band" (see caption Fig. 1) into two subbands between 1.9 to 2.2 /im and 2.2 to 2.5 fim, little 
contrast is apparent in the former while a large contrast is seen in the latter. 

The nature of the state variables chosen for decorrelation has been discussed extensively 
in SVT08-S1. One change is in the array temperature vector: instead of using temperatures 
from the probe at the NIC-1 mounting cup, we use the diode properties of the detector 
pixels as proxy for the detector temperature because the measurable detector bias voltage 
is highly temperature sensitive. These temperatures were extracted by the STScl staff on 
request (courtesy N. Pirzkal, priv. comm.). We find that the peak-to- valley variations of 
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the state variables are similar (i.e. < 0.1 pixel for PSF displacement) to those seen in the 
observations of SVT08 for all orbits included in the fitting procedure, i.e. O2 — O5. A large 
rotation of the dispersion axis is apparent in Oi, and Oi data were discarded in order to 
satisfy our underlying assumption - that spectrophotometry systematics are linear in small 
perturbations of the state variables. After Oi the rotation settles, although O2 shows a larger 
rotation from mean than the final three orbits. 

The framework for decorrelation of the light curves is the joint energy minimization 
procedure described in SVT08, i.e. modelling and removal of the systematics was performed 
simultaneously, with the cclipsc-dcpth, a{\), left as a free variable. The unique minima 
in the a~X space determines the emission spectrum. The point-to-point photometric scatter 
(see Fig. 1), which is stochastic in time but common to all wavelength channels, is estimated 
using optimally weighted fit-residuals as in SVT08-SI. For all channels, the final residual 
vectors are examined via Lomb-Scargle periodograms to ensure that correlated 1// noise, as 
well as noise at the HST orbital period and its overtones, is removed. An alternative scheme, 
that uses just the baseline light curves to fit for a model (SVT08-S1) and then interpolates 
to determine the emission spectrum, gave answers that were statistically consistent with 
the joint energy minimization scheme; this implies that just the baseline orbits provide an 
adequate model. 

The importance of the various state variables in modelling the overall systematics is 
assessed via principal components. The model in each light curve is first decomposed along 
each state vector, which is converted to a channel-wise correlation matrix. Each correlation 
matrix is diagonalized, along with the computation of the principal vectors. We find that a 
maximum of two principal components is required to model the bulk of the variance in the 
channel light curves. The rotation 9, and the displacements x, y, are the most important 
projections onto the first principal component, which models between 55-90 % of the variance 
across channels. The spacecraft orbital phase, (j)H, is the single most important projection 
onto the second principal component. This behavior is largely explained by the fact that 
most elements in the set of state vectors show some dependence on the spacecraft orbital 
period. 

The final emission spectrum (planet/star ratio) is shown in figures 2 & 3 and tabulated 
in Tab. 1. The error bars contain three major components in quadrature. These are as 
follows: (1) the random noise, including detection noise, errors in the estimation of the 
state vectors, and any remaining systematics are contained in the light curve residuals, (2) 
stochastic uncertainties in the model are computed by bootstrapping, and (3) gross errors 
in the model fitting. The latter are estimated either via Monte Carlo methods by applying 
random perturbations to the model components (SVT08-SI), or via the surfaces obtained 
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from the fitting procedure. Finally, we draw attention to the fact that the penultimate bluest 
channel (1.525 /xm) has an unphysical eclipse depth; specifically, the wavelenth channel is 
negative at 2.5cr. This may be due to the fact that this channel contains a dead pixel in 
the brightest part of the spectral image. We have not evaluated the effect of such dead- 
pixels on the spectrophotometry in a strict statistical sense; consequently such effects are 
not reflected in the error-bar for that channel. Previous ground-based ob servations have bee n 
used to estimate the K-band planet/star contrast upper limit as 0.00039 (IBarnes et al 1120071 ). 
We find the contrast in the same passband (averaging the contrast in our spectral channels 
between 1.98 and 2.38 yum) to be 0.0004; given the challenges associated with ground-based 
secondary eclipse measurements, we do not consider the difference significant. 



3. Discussion 



The interpretation of emergent spectra is generally complex because molecular bands 
could appear either in emission, absorption, or both, depending on the detailed tempera- 
ture, pressure, and chemical profile of the planetary atmosphere. Nonetheless, when molec- 
ular species with detectable spectral features are present, radiative transfer models can be 
used to retrieve the temperature and abundance structure of the photosphere of emission. 
This approach, used ex tensively for the study of planetary atmospheres in our solar system 
( iGoody &: Yung Ill989l ). is the basis of our spectral retrieval. 



We undertook a detailed analysis using a line-by-line, radiative-transfer model developed 
for disk- averaged planetary spectra and subsequently adapted for the specific case of hot- 



Jupiters (ITinetti et al. 112006 



2007bl ). The model covers a pressure range from 1 < P < 10~ 



bar in which the individual model layers are populated by abundance profiles for molecules 
used in the retrieval process. The opacity contribution of each molecule is computed based 
on its mixing ratio q, local density (p), and temperature in accordance with the assumed 
T-P profile. Our compendium of line-lists is the same as that given in SVT08. Trial plan- 
etary spectra were converted to con trast spectra u sing an ATLAS-9 model-spectrum of a 
K1-K2V star with solar abundances (IKurucz Ill993l ). To retrieve molecular abundances and 
the temperature profile, we iteratively compared the modelled contrast spectra to the ob - 
served NIR spectra as well as mid-IR Spitzer-IRAC photometry (jCharbonneau et al. II2008I ). 
We found that the IRAC photometry did not provide constraints for molecular mixing ratios, 
although those data were consistent with a decreasing T profile. Thus, the results for molec- 
ular abundances presented here are independent of IRAC data that probe lower pressures in 
the atmosphere. 



The significant spectral modulation present in the near-IR contrast spectrum is best ex- 



- 6 - 



plained by the presence of water vapor and the carbon-bearing molecules CO and CO2. The 
model requires a decreasing temperature with altitude between ~ 1 < P < 10^^ bar (the rel- 
evant photosphere). The modelling began by including the minimum number of molecules 
(H2O, CO, and CH4) composed from the most abundant reactive elements (H,0,and C). 
These molecules are highly plausible; H2O and CH4 have already been inferred via trans- 
mission photometry and spectroscopy (ITinetti et al. Il2007al : ISwain et al. II2008I ). and CO, 
thermochemically the most stabl e carbon molecule on the h ot day-side, has been inferred 
from photometry (IBarman II2008I : ICharbonneau et al. II2008I ). However, this minimal com- 
bination of molecules is unable to reproduce the large contrast deficit in the 1.95-2.15 fim 
region of the spectrum, suggesting additional opacity by the 2 fim combination-bands of 
CO2. The inclusion of CO2 in the model reproduces the major features of the observed spec- 
trum. Fig. 2 illustrates the contribution of individual molecules to the model spectrum; the 
significance of CO2 is evident. This suggests searching for evidence of other CO2 bands, such 



as th e 15 fim band, which may be present in the Spitzer mid-IR spectra (IGrillmair et al. 
2OO7I ). 



Although CO2 is not the primary carrier of carbon in the pressure and temperature 
range of relevance here, it is thermochemically predicted to be fairly abundant (10~ ^ for 



solar abundance) when CO is the major carbon-bearing gas (iLodders fc Feglev 



has been discussed in the co ntext of low-gravity brown dwarf spectral models (ISaumon et al 



2002) and 



2OO3I ; iMainzer et al. 1120071 ). It is formed via the net thermochemical reaction, 



CO + H2O ^ CO2 + H2, 
Additionally, CO2 is readily produced by the reaction 

CO + 0H^C02 + H. 



where the hydroxyl radical OH is derived from photolysis of H2O (ILiang et al. II2003I ). 



While the observed spectrum places strong constraints on what major molecules are 
present, the model fitting must deal with a degeneracy between dT/dP and the molecular 
mixing ratio. In order to estimate the range of plausible mixi ng ratios for the inferred species 



our f inal family of models explored a range of T profiles (IBarman I l2008l : iBurrows et al. 



20081 ). with constant vertical mixing ratios for H2O, CO and CO2, and a vertical mixing 
profile for CH4, as derived by Liang et al. 2003. Based on this family of models, we find the 
following range for mixing ratios at pressure altitudes relevant here: H2O (c ~ 0.1 — 1 x 10~^), 
CO (c~ 1-3x10-^), CO2 (c~ 0.1-1x10-6), andCH4 (c < 1x10"^). These values imply 
a 0.5 < C/0 < 1 together with a metallicity that is potentially subsolar. The inclusion of 
CO2 in our models does not explain the contrast deficits at the band-edges, particularly the 
blue edge at 1.6 yum (Fig. 3). One possibility is the presence in the atmosphere of trace 



- 7- 



amounts of hydrocarbons such as acetylene (C2H2), ethane (C2H6), and possibly ammonia 
(NH3). For example, C2H2, CH4 and C^ Hr can be produce d via photolysis of CO (the 
primary reservoir for carbon) and water (ILiang et al. II2004I ). while ammonia is likely to 
be found in low concentrations (c ^ 10"'') in thermal equilibrium with N2, which is the 
dominant reservoir for nitrogen in these conditions. Indeed, including one or more of the 
above molecules does improve the model fit at the band-edges. However, a firm assertion of 
their presence amongst a myriad of possibilities would require better data and/or broader 
wavelength coverage. 



4. Conclusions 

In summary, we have presented the first near infrared spectrum of light emitted by a 
hot- Jupiter type planet. Using a radiative transfer model, we determine that the molecules 
H2O, CO and CO2 are likely present on the dayside of HD 189733b, and we are able to 
estimate abundances for these species. Although we cannot tightly constrain the slope of 
the temperature profile, the observed absorption features indicate that temperature decreases 
with altitude at pressures between 1 and 0.01 bar. In addition, there are residual features 
in the spectrum that could be explained by including opacities of other species that are 
presently not well constrained. 

This work represents an initial step in exploiting the power of spectral analysis (even at 
low resolutions of i? ~ 40) in determining the chemical composition of extrasolar planetary 
atmospheres; it also illustrates the extraordinary potential of the NICMOS instrument for 
characterizing bright, transiting exoplanets. In a previous paper, we presented a transmis- 
sion spectrum of this planet at the terminator (SVT08), in which methane is seen to be 
more abundant (c = 5 x 10~^) in the higher, cooler regions of the terminator region atmo- 
sphere. However, it is difficult to compare directly the previous terminator results with our 
current dayside results because they probe atmospheric regions with significantly different 
temperatures and altitudes on this highly irradiated planet. The development of sophisti- 
cated global circulation and chemistry models could significantly advance the state of the 
art in this respect. 

We thank Tommy Wiklind, Nor Prizkal, and other members of the Space Telescope 
Science Institute staff for extensive assistance in planning the observations and for providing 
advice about ways in which the observations could be optimized. We also thank Jonathan 
Fortney for helpful recommendations on improving the presentation of this material. G. 
Tinetti was supported by the UK Sciences & Technology Facilities Council. A portion of the 
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Fig. 1. — This figure shows a set of normahzed hght curves plotted against time (in minutes; 
to corresponds to MJD 54220.1442888) in a top-to-bottom sequence. To aid the eye, we 
plot the data together with a set of model light curves calculated assuming a uniform-disk 
for HD 189733b passing behind the parent star. (A) The top curve is a raw broadband 
or "white-light" light curve showing all five orbits. The small broadband eclipse depth of 
2.9 X 10~^ is barely obvious in this light curve. (B) Next is a raw "red-broadband" light 
curve covering a wavelength range of 1.9-2.44 fim. The somewhat larger photometric depth 
of 4.1 X 10~^ is more easily discernible here. Orbit 1 data are poorly behaved over this 
restricted band. (C) Next is the "red-broadband" light curve after the removal of point-to- 
point photometric correlations. This step tightens the scatter in the light curve. Note that 
temporal correlations are not removed and are clearly seen in the light curve. (D) At the 
bottom is a fully processed light curve with all correlations modelled and removed. Gaps 
in the time series occur because HD 189733 is not in the continuous viewing zone for the 
Hubble Space Telescope. 
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Wavelength (micron) 

Fig. 2. — A comparison of the observed spectrum (black markers with error bars) to a 
sequence of component models showing the affect of water together with individual carbon- 
bearing molecules included in our final model. This comparison shows the portion of the 
spectrum where each molecule has significant opacity. In particular, the role of CO2 in 
producing the observed minima between 1.95 and 2.15 /xm is evident. 
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1.6 1.8 2.0 2.2 2.4 

Wavelength (micron) 

Fig. 3. — The near-IR dayside emergent spectrum used in our analysis with ±1 — 0" errors 
shown (black), together with a model spectrum (orange) containing the molecules H2O, CO, 
CH4 and CO2, which are responsible for the absorption features (the strongest of which are 
identified above). The model also includes H2— H2, which contributes continuum opacity. 
The fit residuals suggest that one or more additional molecular species may be present. 
Although the fit is improved slightly by including C2H2, C2H6, or NH3, additional data is 
required to make a strong case for the presence of additional molecular species. 
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Planet/Star Contrast Spectrum 



wavelength 


Fplanet/ Fstar 


error 


2.446 


0.000491 


0.00009 


2.389 


0.000633 


0.00009 


2.332 


0.000592 


0.00008 


2.274 


0.000860 


0.00008 


2.216 


0.000714 


0.00008 


2.159 


0.000452 


0.00007 


2.101 


0.000250 


0.00007 


2.044 


0.000122 


0.00006 


1.986 


0.000030 


0.00007 


1.929 


0.000128 


0.00007 


1.869 


0.000389 


0.00006 


1.814 


0.000331 


0.00006 


1.755 


0.000532 


0.00006 


1.698 


0.000452 


0.00007 


1.640 


0.000177 


0.00007 


1.584 


0.000021 


0.00006 


1.525 


-0.000190 


0.00007 


1.476 


-0.000089 


0.00008 



Table 1: Data and l-o" error for the contrast Fpianet/ Fstar spectrum of the dayside of 
HD 189733b. 



